Discrete thinning dynamics in a continuum model of metallic nanowires 
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Transmission electron microscopy experiments have recently observed gold metal nanocylinders 
to thin down in a discrete manner: A kink — a step of order of one atomic layer — nucleates at one 
end and then moves along the wire, leaving a thinner cylinder behind it. In this paper, I show that 
a similar thinning process takes place within the nanoscale free-electron model, a structural and 
dynamical model of nanowires that treats the electron-confinement effects exactly while replacing 
the atomic structure by a continuum. Electron-shell effects, previously shown to be responsible for 
the stability of wires with magic radii, favor the formation of kinks connecting magic cylinders. A rich 
kink dynamics including interkink interactions ensues and is similar to that observed experimentally. 

PACS numbers: 68.65.La, 66.30.Pa 
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The formation and thinning dynamics of metal- 
lic nanowires suspended between macroscopic con- 
tacts has recently been imaged, in real time and 
with atomic resolution, using transmission electron mi- 
croscopy {TYM)X^^^^^ Almost perfect Au (Refs. H, i, 
and d) and Ag (Refs. d, and 0) cylindrical wires, with 
diameters ranging^ from 5 to 15 A, are found to form 
and be stable for seconds at room temperature. These 
cylinders usually do not break up at once but are instead 
seen to thin down one step at a time by an amount close 
to one atomic layer.— Real-time movies^ of the thinning 
dynamics reveal that wires can remain in a metastable 
state for a time of order of a second, until a kinklike struc- 
ture nucleates at one end of the wire and subsequently 
moves along the cylinder. The kink formation and initial 
displacement over a significant portion of the wire are ex- 
tremely rapid, of the order of the time resolution of the 
experiment, while subsequent motion is slower and takes 
a significant fraction of a second^ for a displacement of 
order of 1 nm. In some cases, a kink appears to stop 
along the wire until it is joined by another one, at which 
point they move along together until they are absorbed 
by the contact at the other end of the wire. 

Due to the large fraction of surface atoms — with low 
coordination numbers — in such nanowires, surface effects 
are particularly important: Barring an additional stabi- 
lizing mechanism, they are expected to trigger a Rayleigh 
instabilityii^ that would break the wire up. Surface ten- 
sion indeed seems a likely candidate for the force driving 
nanowire thinning but cannot account for the step by 
step nature of the process. A layer by layer thinning of a 
crystalline structure could, in principle, account for itf^ 
but comparison of TEM images with simulations^' leads 
to the conclusion that such thin nanowires do not have 
a crystalline structure. They adopt instead a multishell 
structure with cylindrical symmetry, similar to carbon 
nanotubes, which would not be classically stable. This 
suggests that the atoms adjust themselves as closely as 
possible to the electronic structure, which favors cylin- 
drical symmetryfii which cannot be done while keeping 
a crystalline structure. 

Electron-shell effects — similar to those well-known in 



cluster physicsi^ — have been shown to be essential to the 
wire stability^ '^'^1^^ providing an energy barrier sufficient 
to make nanowires metastable, with lifetimes^^ of order 
of a second. This energy barrier arises from the competi- 
tion between surface tension and electron shell effects, the 
latter resulting from the quantum confinement of conduc- 
tion electrons within the cross section of the wire. The 
"jerkiness" of the thinning suggests a thermally activated 
process, where random fluctuations are occasionally large 
enough to overcome the barrier, thus triggering thinning 
to the next metastable wire.— This picture is reinforced 
by the fact that thinning forced by electron irradiation of 
the wire follows the same mechanism, but is much more 
regular. 

Any meaningful theoretical description of this prob- 
lem requires a model that captures both surface and 
electron-shell effects. In addition, it needs to be simple 
enough to allow simulations of large systems over long 
time scales. Both of these conditions are fulfilled by the 
nanoscale free-electron model (NFEM),^^'^^'^* a contin- 
uum model where the atomic structure is replaced by a 
uniform, positively charged background, and the empha- 
sis is put on the electronic structure. Ah initio methods, 
which include both electronic and ionic structures, sat- 
isfy the first of the above conditions, but not the second. 
They have been used to confirm the atomistic shell struc- 
ture of gold nanowires through energy minimization^i^ 
but they are computationally quite intensive and lim- 
ited to relatively thin wires even for structure calcula- 
tions. As the observed thinning dynamics involves hun- 
dreds, if not thousands, of atoms and takes place over 
a sizable fraction of a second, it is far beyond the reach 
of ab initio simulations, ^^'^"^ which are typically limited 
to a few dozen atoms and cover a nanosecond. Even 
classical molecular dynamics j^i^i which do not include 
electron-shell effects, can only simulate a few nanosec- 
onds at most. 

The NFEM has successfully explained the linear sta- 
bility of cylindrical nanowires with "magic" radii as re- 
sulting from a competition between surface and electron- 
shell effects.— As a continuum model, the NFEM can 
obviously not address issues such as the infiuence of the 
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crystalline structure of the macroscopic contacts on the 
wire stability*^ The model has proved useful in under- 
standing ionic dynamics through surface self-diffusion 
over long time scales. "'^^ In particular, a random wire has 
been shown to naturally evolve into a universal equilib- 
rium shape consisting of a perfect cylinder of a magic 
radius, connected to thicker leads through an abrupt 
junction] ^^1^^ A stochastic model of the dynamics of 
a nanocylinder under thermal fluctuationsi^ has pre- 
dicted long lifetimes compatible with experimental obser- 
vations. The escape mechanism for long wires has been 
shown to consist of the nucleation of a kink at one end 
of the wire. The kink subsequently propagates along the 
wire, leaving behind a thinner or thicker cylinder. This 
suggests that a dynamics similar to experiments may take 
place within the NFEM. 

Motivated by these results, I argue that, despite be- 
ing a fully continuum model, the NFEM contains an in- 
trinsic length scale — the Fermi wavelength Xp — which, 
through electron-confinement effects, provides some de- 
gree of "discreteness" for the ionic dynamics. The com- 
petition between surface and electron-shell effects favors 
the formation of kinks, or solitons, connecting cylinders 
of magic radii. The dynamics of these solitons is shown 
to be in semiquantitative agreement with experiments in 
many details. In particular, the model predicts an ex- 
tremely fast nucleation and initial propagation of kinks 
at the wire end, followed by a regular, slower motion. 
Interactions between kinks, which can be either attrac- 
tive or repulsive, profoundly affect their dynamics. For 
some wires, a soliton is seen to be essentially stopped by 
its interaction with another one, until they combine and 
move as a single, larger kink, much like what is observed 
exp er iment allyi^ 

The paper is organized as follows. In SecUl the NFEM 
is introduced and some of its main results directly rele- 
vant to the present article are summarized. The forma- 
tion and propagation of a soliton are discussed in Sec. [Til 
while interactions between multiple kinks are analyzed in 
Sec. mil A discussion of the results and their relation to 
existing experiments is provided in the last section. 



I. THE NANOSCALE FREE-ELECTRON 
MODEL 

The NFEM is a continuum model of open metallic 
nanosystems with emphasis on the electronic structure, 
which is treated exactlyJi It is thus particularly suit- 
able as a model of metal nanowires, where electron- 
confinement effects have been shown to influence both 
transport and cohesive properties^^i^ (see Ref. [1^ for 
a recent review.) The discrete atomic structure is re- 
placed by a uniform, positively charged background (jel- 
lium), which provides a confining potential for electrons. 
Electronic degrees of freedom are described using a free- 
electron model, thus neglecting interactions, except inas- 
much as they rescalei^*^*^ macroscopic quantities such 



as the bulk energy density ub, and the surface tension 

This model is particularly suitable for simple metals, 
with good screening and a close-to-spherical Fermi sur- 
face, such as those with a single s-electron conduction 
band at the Fermi surface. Such conditions are fulfilled 
for alkali metals such as sodium and for noble metals such 
as gold and silver, although d electrons may play a role for 
noble metals. The NFEM has provided an understand- 
ing in simple physical terms of many transpor t^^'^^i^^i^^ , 
stabilityj^iiiii^i^ and dynamica l^'^'^^i^^'^^ properties of 
alkali and noble metal nanowires, in quantitative agree- 
ment with experiments.. 

While the NFEM has no such restriction, only axisym- 
metric wires are considered in this paper. This choice is 
justified by the facts that (i) the most stable wires are 
axisymmetric,"^^^ with only a small fraction of low con- 
ductance wires having a broken axial symmetry, and (ii) 
the dynamics tends to decrease surface area and thus fur- 
ther favors axisymmetric wires. In addition, axial sym- 
metry greatly simplifies the numerical treatment of the 
dynamics, as the wire shape can be described by a single 
radius function R{z,t). 

A nanowire connected to macroscopic contacts be- 
ing an open system, the electronic energy is given by 
the grand canonical potential fJe. Like any extensive 
thermodynamic quantity, can be written as a Weyl 
expansion's. — a series in geometrical quantities such as 
system volume V and surface area S — complemented by 
a mesoscopic, fluctuating contribution 

ne[R{z)] = ljbV + asS + sn, (1) 

where the values oi ujb and as may be chosen to match 
the bulk properties of the metal to be described. As 
results are independent of ujb, its free-electron value 
ujb = —iEpkp/lbn'^ is used, while the surface tension 
is set to (Ts = 1.256 N/m, a value appropriate for the 
description of goldi^ 

Assuming the wire cross section varies slowly along the 
wire (adiabatic approximation) , the mesoscopic contribu- 
tion, which is of purely quantum-mechanical origin, may 
be written as 

Sn[R{z)] ^ [ dzVsi,eii[R{z)], (2) 
Jo 

where the electron-shell potential Vsheii{R) can be com- 
puted using a semi-classical approximation^i^ The name 
"electron-shell potential" refers to the analogous shell 
closing observed in metal clusters. ""^^ Although there is 
strictly speaking no shell closing in metal nanowiresji^ 
since they are open systems, a large gap in the transverse 
eigenenergies at the Fermi energy increases the stability^S 
of the corresponding nanowire. This results in a deep 
minimum of Vsheii (R) j depicted in Fig. [T]Ja) for a cylin- 
drical wire as a function of its radius R. A linear stabil- 
ity analysis2ii2ii22 shows that wires within a finite interval 
around the magic radii are stable toward all small per- 
turbations. 
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FIG. 1: (a) Electron-shell potential VsheiiiR), and (b) atomic 
chemical potential ^cyi{R) for a cylindrical wire, as a func- 
tion of its dimensionless radius kpR, where kp is the Fermi 
wavevector. (The radius range is limited to that correspond- 
ing to the simulations presented here. For a graph over a 
more extended range, see Ref. Q-) The top axis shows the 
conductance values of linearly stable cylinders in units of the 
conductance quantum. Go = 2e^ /h. 



The ionic dynamics is taken to be classical and can be 
assumed to occur mainly through surface self-diffusion, 
as most atoms in thin metal wires are surface atoms j^i^ 
The evolution equation for the radius function R{z, t) 
derives from ionic mass conservation 



Va dt 



+ ■T^J,{z,t) = 0, 

oz 



(3) 



where Va — Sn'^/kp is the volume of an atom in a mono- 
valent metal, and the z component of the total surface 
current is given by Fick's law: 



psDs 27tR{z, t) 
' keT ^1 + [d,RY dz ' 



(4) 



Here, ps and Ds are, respectively, the surface density of 
ions and the surface self-diffusion coefficient, and dzR = 
dR/dz. 

The ionic chemical potential fi[R{z)] can be computed 
from the energy change due to the local addition of 
the volume Va of an atom to the system. Within the 
Born-Oppenheimer approximation, and assuming that 
the electrons act as an incompressible fluidfi^ii^i^ the 
chemical potential is given^* by the functional derivative 
fi[Riz)] = [Va/(27ri?)] X [Srie/SRiz)]. Starting from Eqs. 
(P) and one obtains 



fi[R{z)] = ^0 



Va / 2<jsdC[R{z)] ^ dVs 



shell 



2ttR 



v/1 + idzRY 



OR 



(5) 



where po 
dC[R{z)] = TT 



uJsVa is the bulk chemical potential. Here, 



i +(d^R)'^ I is the local mean curvature 
of the wire and results from the functional derivative of 



the surface term. The chemical potential of a cylinder 
l^cyi{R) = l^[R{z) = R] is plotted as a function of the 
radius R in Fig. [Ijb). 

The precise value of Ds in Eq. ([4]) is not known for 
most metals, but it can be removed from the evolution 
equation by rescaling time to the dimensionless variable 
T = ujotj with the characteristic temperature-dependent 
frequency luq = psDsTp/T. For comparison with exper- 
imental time scales, one can estimate that for quasi-one- 
dimcnsional diffusion, Dg ~ i^ua^ exp(— iSs/fc^r), where 
i^D is the Dcbyc frequency, a is the lattice spacing, and 
Es is an activation energy comparable to the energy of a 
single bond in the solid. 

As mentioned above, a linear stability analysis® finds 
intervals of stable radii. The criterion for linear stability 
of a cylindrical wire has been showni^ to be 



dpcyi (R) 
dR 



> 0. 



(6) 



A linearized dynamical theory^^ suggests that unstable 
wires develop an exponentially growing instability with a 
well-defined wavelength corresponding to the maximally 
unstable mode. This instability was argued to saturate 
and eventually lead to a phase separation of the wire 
into thick and thin segments of stable radii. Simulations 
using the full dynamics, defined by Eqs. ([3]) -([5]), have 
confirmed this and have shown that the phase separation 
occurs via a complex dynamics involving kink interac- 
tions and annihilation^^ 

Note that, despite the apparent simplicity of the evolu- 
tion equations [Eqs. ((21)^®] as written above, the result- 
ing partial differential equation for t) is fourth order 
in z derivatives and highly nonlinear. Its classical coun- 
terpart, with Vsheii = 0, has been extensively studied.— 
It was shown to have stationary states corresponding to 
shapes of constant mean curvature, the sphere, the cylin- 
der, and the unduloid of revolution, the latter being al- 
ways unstable. The addition of the electron-shell po- 
tential, of quantum-mechanical origin, stabilizes magic 
cylinders connected to unduloid-shaped leads as the uni- 
versal equilibrium shap c"*! ^'— ' and leads to the rich kink 
dynamics, as discussed in this paper. 

Given the universality of the equilibrium shape and 
the "sharpness" of the connection between the unduloid- 
shaped leads and the cylindrical wiie^^^ one can focus 
on the dynamics of the cylindrical part of the wire, as- 
suming that the lead is large enough to serve as a source 
or sink of atoms for the wire without undergoing signif- 
icant changes. In this case, Neumann boundary condi- 
tions, dzR — at both ends of the wire, z — and 
z — L, are found to best describe the connection to the 
leads. The diffusive dynamics being volume conserving, 
a sink or source of atoms at one or both wire ends must 
be explicitly added to mimic a nonequilibrium situation, 
where an ionic current flows between the cylinder and the 
leads, thus allowing for thinning or growth of the wire. 
This is done by adding a boundary current term to the 



right-hand side of Eq. Q : 

Jbnd = JlS{z) - JrS{z ~ L), (7) 

where a positive J;(r) is a source of atoms at the left 
(right) wire end, while negative values correspond to 
sinks of atoms. 

II. SINGLE KINK DYNAMICS 

Motivated by the experimental observation of wire 
thinningj^ I will consider mainly the case of a single sink 
of atoms at the left end of the wire, J; < and Jr — 0. 
A number of simulations for various wire lengths and 
boundary currents have been performed. Figures [21 SI 
and O show results for kpL ~ 150, or L ~ 12 nm for Au, 
with a boundary current J; = — IO^^ttwo. The parame- 
ter dependence of the results, as well as the influence of 
a second sink of atoms, or fluctuations of the boundary 
current, are discussed in Sec. IIIII 

A boundary current J; < creates, following Eq. ([4]), 
a positive gradient of chemical potential that pro- 
gressively expands along the wire. From Eq. ([5]), a stable 
wire corresponds to a positive slope of ^cyi 

(R). Thus, 

a positive gradient of R{z) corresponds to the gradient 
of Once the radius change at the boundary is 

large enough to cross a stability threshold, correspond- 
ing to a minimum of ^cyi{R), an instability sets in^'^'^* 
and leads to the formation of a soliton that will prop- 
agate along the wire. In most cases, like the transition 
from G/Go = 42 ^ 34 (Fig. d]), or G/Gq = 23 ^ 17 
(Go — jh is the quantum of conductance), a single 
kink nucleates and propagates along the whole wire. The 
present section discusses this case in detail. Sometimes, 
however, a second kink nucleates before the first one has 
reached the end of the wire, as is the case in the transi- 
tion G/Go = 34 ^ 27 ^ 23 (Fig.H) If this happens, the 
two solitons interact attractively and eventually combine 
and propagate as a single, larger kink (Figs. [Hand [5]) 
This is the simplest case of kink interactions, which are 
discussed in Sec. IIIII 

FigurelUa) displays the radius function R{z,t) at dif- 
ferent evolution times (with alternating solid and dashed 
lines for clarity), with a dotted horizontal line showing 
the channel opening threshold, used to track the kink 
position Zk- The radii corresponding to unstable cylin- 
ders are emphasized with a gray background. The corre- 
sponding chemical potential profiles ii{z,t) for the first 
four configurations are shown in Fig.[2{b). (The last two 
configurations are omitted as they are indistinguishable 
from the last displayed one, which looks constant on the 
scale used.) The lower panel (c) shows function 
of time (solid line), both on linear (main plot) and log- 
log (inset) scales, as well as the kink velocity Vk (dashed 
line). One can clearly identify two phases with different 
time evolutions of the kink position: the kink formation, 
during which Zk increases quickly but the kink velocity 
Vk decreases abruptly, and the kink propagation, when 
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FIG. 2: (a) Radius function R{z) at various evolution times 
during the transition G/Go =42 -^34 for a wire subject to 
a sink of atoms at 2 = 0. Solid and dashed lines are used on 
alternate curves for clarity, while the gray areas correspond 
to unstable radii, and the dotted line shows the radius used 
to define the kink position Zk- (b) Corresponding chemical 
potential profiles fi{z) for the first four curves in (a). The 
last one looks constant on this scale, but has a slight positive 
slope left of the kink, (c) Kink position Zk (solid line) as a 
function of time on linear and log-log (inset) scales. The kink 
position for the shapes plotted in the top panel are marked 
with corresponding symbols on each graph. The dashed line 
gives the kink velocity Vk, with units given on the right axis. 



the kink moves without deformation at a constant veloc- 
ity vq. The kink formation actually takes place in two 
stages, as will be discussed shortly. 

As soon as the wire radius at the sink of atoms reaches 
a critical radius i?c, corresponding to a marginally sta- 
ble wire [upper limit of gray area in Fig. [H^a)], for which 
d^cyi I dR{Rc) = (first curve from the left), an instabil- 
ity sets in^ » and grows exponentially, leading to the quick 
formation of a kink (second curve from the left). During 
this first stage of the kink formation, the radius "over- 
shoots" ahead of the kink and creates a "bump" visible 
in the first dashed curve. In a second stage, this over- 
shoot progressively spreads over the whole wire (second 
to fourth curves from the left), dragging the kink behind 
itself at a quickly reducing speed. The kink subsequently 
propagates without deformation (remaining curves) with 
a constant velocity, proportional to the boundary current 
I J; I , provided the wire is stable enough, not too long, and 
the current small enough (see Sec. IIII[) . 
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A. Kink formation 

When the wire radius at the boundary reaches an in- 
stabihty threshold i?c, an instabihty with a single wave- 
length A = 27r/q starts growing exponentially at a rate 
Wfi^ as discussed in Sec. HI As the instability is driven 
by the boundary current, it is initially localized at and 
propagates from the wire end. This is accounted for by a 
propagating envelope function chosen as an exponential 
with decay length ^ and propagation speed v, which mul- 
tiplies the perturbation function SR{z), which can thus 
be written as 

SR - Re {exp[i{ujt - qz)] exp[-(z - vt)/^]} . (8) 

Although there is no well-defined kink at this point, its 
position can be defined as the point where SR reaches 
a fixed value i5. Results depend weakly on that valuer^ 
which is chosen to be the threshold for conduction-band 
closing. 

As happens for the instability of a cylinder j^'^i^^ the 
growth saturates, leaving a well-formed kink, once it 
reaches a zone of stable radii, i.e., grows out of the gray 
areas in Fig. [2^a). In this particular case, the instability 
grows through two unstable zones, shown by the two gray 
areas, as the intermediate radii are not stable enough to 
stop its exponential growth. 

The product of the decay length ^ of the envelope and 
wave vector q of the perturbation can be shown^^l using 
the linearized version of evolution equations [Eqs. ([3])- 
([5])], to be a constant ~ 3.2, so that the instability is 
only visible over essentially one wavelength A. Its first 
half wavelength becomes the kink, while the second half 
wavelength, much reduced in amplitude, is responsible 
for the observed overshoot. 

Solving SR{zk,t) = S for small qzk{t), neglecting the 
envelope at first, the kink position during its formation 
is found to be 

kpzkir) (X \/ujo{t - ^o), (9) 

where to is the kink formation time. This result is in- 
dependent of the boundary current, as long as the latter 
is sufficiently small compared to the growth rate of the 
instability. Corrections to this behavior due to the enve- 
lope are small in its range of validity, as its propagation 
speed is much smaller than its growth rate. This t^^"^ be- 
havior of the kink position is observed only in the early 
stages of the kink formation, At — ujo{t — io) ^ 2 x 10"*, 
as can be seen in the inset of Fig. [21 where has been 
plotted as a dotted line for comparison. 

B. Overshoot spreading 

The progressive spreading of the overshoot over the 
whole wire, driven by its far-from-equilibrium chemical 
potential [rightmost solid line in Fig. [DJb)], speeds up 
the kink propagation, as is discussed in this section. 



(a) (b) 




FIG. 3: Schematic diagram of the radius function R{z), and 
corresponding chemical potential ^[z), during (a) the over- 
shoot propagation ahead of a kink, (b) a soliton propagation, 
(c) interaction of two kinks, and (d) interaction of a kink and 
antikink. The diagrams define the radii and lengths used in 
the text, as well as the chemical potential drop A/i. The 
arrows labeled Jk mark the positions of sinks or sources of 
atoms. In panel (d), only the radius function is shown, and 
the dashed line represents R{z) before nucleation of the up- 
ward kink. 



Prior to the kink formation, the wire radius is driven 
by the boundary current out of its most stable value, 
which corresponds to a minimum of Vsheii (R) ■ It is in- 
stead close to its marginally stable value which, through 
Eq. ([S]), corresponds to a minimum of fj,cyi{R). The over- 
shoot brings a portion of the wire back to its most stable 
radius Rr, creating a negative gradient of chemical po- 
tential and, according to Eq. a corresponding ionic 
current in the forward direction. This current favors the 
spreading of the overshoot and drags the kink along with 
it. 

The chemical potential along the wire during overshoot 
propagation is depicted schematically in Fig.[2{a). It has 
a constant slope, corresponding to a slight gradient of 
R{z) (see beginning of the present section), proportional 
to the boundary current on the left of the kink and is 
constant in the cylindrical part of the wire on the right 
of the overshoot. It is essentially constant over the kink 
and drops by an amount A/z < over the extent of the 
overshoot. As discussed above, A/j, is essentially fixed 
by the potential difference between the most stable and 
the corresponding marginally stable wires in the stabil- 
ity interval considered, and is thus constant during evo- 
lution. Consequently, the gradient of potential A/z/Lq 
over the overshoot, of length Lq, decreases (in absolute 
value) with time. 

With these elements in mind, the overshoot spread- 
ing, and the corresponding kink propagation, can be 
fully understood from the fact that the evolution [Eq. 
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Q] is volume conserving. Let us consider a portion of 
the wire of length Lg > Lq, starting at the overshoot 
maximum. Approximating the overshoot as a cone of 
average radius Rr — AR/2 [see Fig. EJa)], the volume 
Vo = TT[{Rr - AR/2)'^Lo + {Rr - ARf{Lo - io)] of this 
section of the wire increases over time due to an incoming 
current J = — 27rwo^rA/i/Lo on its left-hand side, while 
there is no current on its right-hand side, which forces 
the overshoot expansion. Integrating Eq. ([3]) along this 
portion of the wire, one obtains the following mass con- 
servation equation: 



Va dt ^ -^'^f' 



J, 



right 



J, 



(10) 



where J; 



left(right) 



are the ionic currents on the left- 



hand (right-hand) side of the section of the wire consid- 
ered; in this case, their values are J and 0, respectively. 
This equation yields the ti me dependence of the over- 
shoot length, Lo(At) — loV At, where At measures the 
time from its formation, and 



AVaRr^H 



EpARiRr - 3AR/4) ' 



(11) 



The parameters in Eq. pT|) . Rr, AR, and Afi — 
IJ^cyi{.Ri) — IJ-cyi{Rr — Ai?) , Can all be obtained from the 
shape R(z,t) at any given time during overshoot propa- 
gation and determine its expansion. 

The time evolution of the kink displacement Azfc(AT — 
T—Ti) = Zfc(T) — z/c(ti), Zk being the kink position, during 
the overshoot propagation, derives from mass conserva- 
tion [Eq. (fTU|) ] for the whole wire and is 



Azfe(AT) 



ARiRr-'-^%Ary^-y^Ar 



{Rr - ARy 



(12) 



In the present simulation, Ji/llIo ^ 1, so that the cor- 
responding term in Eq. p2p may be neglected and only 
the V At behavior is observable, as can be seen in the 
inset of Fig. He) for 10^ < t < 10^ 

How far the kink moves during the overshoot expan- 
sion, as well as how long it takes for the overshoot to 
reach the wire end, can be estimated from Azfe(ATo) -|- 
Lo{Ato) ^ L. It is found that Azfc(ATo) ~ r/L, with a 
proportionality constant 77 depending on the coefficients 
of At and \/At in Eq. For most kinks, 77 0.1-0.3. 
The propagation time Atq is also proportional to L and 
is found to be of order lO'^ — lO''. For the kink of Fig. [21 
the constant is ?7 ~ 0.14, so that the kink jumps ahead 
by about kpAL ~ 20 during this phase that lasts a time 
Ato ~ 4 • 10^, in agreement with the simulation. 



C. Kink propagation 

Once the overshoot has reached the end of the wire, the 
chemical potential profile becomes simpler and is then 
constant for z > Zk{T), as depicted schematically in Fig. 



|3l^b). The same volume conservation argument, applied 
to the whole wire, shows that the kink moves with a 
constant speed determined by the radii on both ends of 
the wire and proportional to the boundary current: 



Zfc(T) = Zfc(Tl) - 



Va 



Jl 



7r(i?2 - Rf) ujo 



{r-r,) 



(13) 



where ti is the time at which the overshoot has been fully 
absorbed at the wire end. The kink velocity obtained 
from Eq. (|13p is in good agreement with that of Fig. 

He). 

None of this discussion depends critically on the sign 
of J;, and the same kind of soliton formation and propa- 
gation is observed for a source of atoms. 



III. KINK INTERACTIONS 

The likelihood of a second kink being nucleated before 
the first reaches the wire end increases with the boundary 
current | J/ 1 and wire length L but decreases with the wire 
stability. As mentioned previously, the chemical poten- 
tial profile on the left of the kink has a constant gradient, 
proportional to J;, corresponding to a small gradient of 
R{z). Inevitably, if the kink were propagating into an 
infinite wire, there would be a time t„ when the kink 
position Zk{Tn) is such that the potential difference due 
to this gradient is large enough for a new instability to 
set in, and thus for a new kink to nucleate. If the wire 
length L is larger than Lc = Zk{Tn), this will happen even 
for a finite wire. The critical length Lc increases with the 
wire stability: A given drop A^ of the chemical potential 
may be enough to trigger nucleation of a new kink for a 
given wire, say, one with G = 27 Gq, but not for a more 
stable wire, such as one with G = 34G'o (see Fig.[T]). As 
the gradient of chemical potential is proportional to the 
boundary current, according to Eq. one clearly has 
Lc oc I J;|. 

During the transition from G/Gq ~ 34 to 23, via 27, 
Lc happens to be shorter than the wire length in the 
simulation discussed above, so that a second soliton is 
nucleated. Figure [3] traces the positions of the two inter- 
acting kinks. Under the influence of the newly nucleated 
kink, appearing around t = 10^ (dashed line), the exist- 
ing kink slows down and even reverses its motion (solid 
line). The approach speed of the two kinks is observed to 
increase as they get closer (see inset). After they inter- 
act, both solitons appear to move in parallel at a constant 
speed. In fact, following the radius function R{z, t) dur- 
ing evolution, plotted vertically with a rightward shift 
proportional to t in Fig. [5l reveals that the two kinks 
combine to form a larger one. 

The interaction between any two kinks can be under- 
stood in simple terms when both are well formed and 
separated enough, i.e., each of them propagates without 
deformation and all overshoot expansion is over. In that 
case, the kinks connect approximately cylindrical seg- 
ments of the wire (each with a slight slope corresponding 
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to a small ionic current) with respective radii of Ri , Rk , 
and Rr, as depicted schematically in Fig.[3]^c). According 
to Eq. (HI, the chemical potential profile on either side of 
the pair of solitons has slopes proportional to the currents 
Jo and Ji on the left- and right-hand sides, respectively, 
which are assumed to be known. There is, in general, a 
linear change of the potential A/i = l^-cyiiRi) — ^J■cyl{Rr) 
over the segment of the wire of length D joining the two 
kinks. Considering mass conservation [Eq. (fTO)) ] sepa- 
rately for two portions of wire containing, respectively, 
only the first and second solitons, and combining the two 
equations thus obtained, the following differential equa- 
tion for the distance D between the kinks is obtained: 

TT dD{t) _ Jo Ji 

27rc^o Rk{Rj-Rf) A^, 
^ {R^,^Rl){Rl-Rf)D{t)- ^ ^ 

As long as the outside currents Jq and Ji are small 
enough, the last term in Eq. dominates, so that 

D^{t) (X {t — to). The proportionality constant can have 
either sign depending on the balance of chemical poten- 
tials of the outside cylinders, as well as the type of both 
kinks (upward or downward solitons). The interaction 
force between kinks can thus be attractive or repulsive. 
It can be calculated as = — ^jf |y, where Vie, given by 
Eq. (P), is the energy of a portion of the wire containing 
both kinks. The result is 
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FIG. 4: Evolution of the kink positions for two interacting 
kinks in a wire subjected to a sink of atoms at its left bound- 
ary. The wire radius before nucleation of the first kink (r < 0) 
is kpR = 12.75, corresponding to a conductance G = 34Go, 
and the wire thins down to kpR ~ 10.7, and G = 23Go, by 
the time the 'combined' kink has reached the wire's end. The 
inset zooms in on the time region where the two kinks interact 
(r 10«). 
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FIG. 5: Radius function R{z) at equidistant times r, show- 
ing the nucleation of a second kink (dashed line in Fig. 3)) 
that interacts and fuses with the preexisting kink around 
r ~ 1.016 • 10®. R{z) is plotted vertically, with a rightward 
shift proportional to the evolution time, for better comparison 
with the inset of Fig. |4l 



where e{R) — 2Tr(JsR + Vsheii{R) is the linear energy den- 
sity for a cylinder of radius R. The interkink force is thus 
constant, independent of their separation. 

Coming back to the solitons of Fig. Hj their interaction 
is found to be attractive, as expected from their observed 
behavior. Even for the largest kink separation, of order 
kpD ^ 100, the last term of Eq. (fT4|l is still one order 
of magnitude larger than the other terms, which may 
thus be neglected. After solving for D{t), the time evo- 
lution for the kink positions z^^ (r) and z^^ (t) can be ex- 
tracted by considering volume conservation for the whole 
wire. Results obtained through this simplified evolution 
are found to be in agreement with the full simulation. 

A simple case of repulsive interaction between two 
kinks is obtained by replacing the sink of atoms by 
a source, J; > 0. While the slightly stable wire at 
G — 30 Go was skipped during wire thinning, or more ex- 
actly a new kink was nucleated so quickly that that radius 
was only observable for a very short time, it is observed 
during wire growth. When the kink to G = 34 Go is nu- 
cleated, its amplitude is large enough that the overshoot 
ahead of the kink — which is now toward lower radii — is 
sufficient to trigger a new instability. This forms an ad- 
ditional kink ahead of the existing one, as schematically 
shown in Fig. [3jd). In that particular case, the dynam- 
ics is found to be dominated by the interkink force [Eq. 
(jlSp ] only in the very early stages of evolution, while 
the external current quickly provides deviations from the 
t^/^ behavior of the interkink distance. Interestingly, the 
contribution of the outside current Jo = J; is attractive, 
so that the sign of the global force may change when 
the kinks are far apart. As the system is overdamped, 
the interkink distance actually relaxes to an equilibrium 
distance, where the force vanishes (as long as no other 
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soliton is nucleated.) 

The general evolution of a system with multiple kinks 
and/or antikinks can be predicted by considering each 
pair of neighboring solitons separately. A set of dif- 
ferential equations similar to Eq. (fT4|) is thus obtained 
for each interkink distance and may be solved, together 
with an additional equation from global volume conser- 
vation, for each soliton position. A general trend of the 
kink dynamics is to decrease their number and increase 
the length of cylindrical segments, as has been observed 
numericall y^^'^^ in the equilibration of a random wire, 
leading to the universal equilibrium shape and the long- 
time evolution of an unstable wire. In the latter case, 
the wire radius switches back and forth between two 
neighboring stable radii, so that Eq. (fTS]) simplifies to 
F = e{Ro) — s{Rk), where Ro and Rk are, respectively, 
the radii outside and between the two kinks considered. 
When F < 0, the cylinder between the two kinks behaves 
exactly as a false vacuum. 



DISCUSSION AND CONCLUSIONS 

The kink dynamics described in Sec. |TT] and IIIII is ac- 
cessible to experiments and has, in fact, already been 
observed for gold^'^^ and possibly for silver— wires in 
TEM experiments: Movies of the real-time dynamics of 
gold wires'' show a stepwise thinning of a cylindrical wire 
from seven down to four atomic diameters through mo- 
tion of kinks along the wire. Solitons are seen to nucleate 
and move very rapidly — faster than the time resolution of 
the experiment — across a significant portion of the wire, 
in agreement with the fast initial propagation predicted 
during kink formation and overshoot expansion in the 
NFEM. Some kinks are seen to stop along the way, as 
during the transition from six to four atomic diameters. 
A second kink appears to join the first one before they 
move along. All this is in semiquantitative agreement 
with simulations presented here as the starting wire with 
a conductance G = 42Go (Fig. ^ has a diameter close 
to seven atomic diameters, while a wire with a conduc- 
tance of 17Go corresponds to five atomic diameters. Note 
that conductance values quoted here are quantized val- 
ues corresponding to ideal, balhstic transport. Measured 
conductances^ are likely to be lower due to backscatter- 
ing of electrons^ from wire imperfections and disorder 
in the leads. 

The simulation includes a sink of atoms, modeled as 
a constant boundary current. This is meant to mimic 
a nonequilibrium situation, where the connection to the 
macroscopic leads allows for atoms to diffuse back and 
forth between the wire and the leads. Experimentally, 
the wires are observed to thin down, presumably due to 
thermal fluctuations, electron irradiation from the elec- 
tron microscope, and/or tension on the wire, thus the 
choice of a sink rather than a source of atoms. A single 
sink of atoms is used in order to simplify the dynamics. 
When two sinks of atoms are included, the dynamics is 



found to be dominated by the larger one: Kinks nucle- 
ate at the position of the larger boundary current and 
propagate through most of the wire before the other sink 
nucleates the corresponding antikink. The thinning dy- 
namics with two unequal sinks of atoms is thus essentially 
equivalent to that of a single sink, except for the absorp- 
tion at the wire end. Finally, thermal fluctuations of 
the boundary current would provide fluctuations around 
the ideal dynamics discussed here but are not expected 
to alter it significantly. The effect of thermal fiuctua- 
tions of the whole wire structure has been studied using 
a stochastic dynamical model.- - The escape mechanism 
from long, metastable cylinders is found to be the nucle- 
ation of a kink at the wire end that propagates along the 
cylinder. The effect of thermal fluctuations may thus be 
included to lowest order as a random fluctuation of the 
boundary current, which would make the thinning more 
"jerky" and thus even more similar to experiments.^'^ 

Material dependence is included in the model by ad- 
justing the surface tension ag in Eq. ([T|) to the appro- 
priate bulk value liiii^ii^ The diffusive dynamics depends 
weakly on this parameter, except for a decrease of the 
critical length Lc for nucleation of new solitons (see Sec. 
HIT)) with increasing surface tension. 

The NFEM is a continuum model that exhibits some 
degree of discreteness, apparent through the formation of 
kinks connecting cylindrical parts of the wire. Although 
the dynamics of the ionic background is fully classical in 
the NFEM, it is strongly influenced by the electron-shell 
potential, which is a result of the quantum confinement 
of the transverse motion of the electrons within the wire. 
The latter provides a "quantization" of the continuous 
ionic structure on a scale given by the Fermi wavelength 
Xp and is responsible for the rich kink dynamics. 

The NFEM assumes that the atomic structure of the 
wire adapts itself to the shape dictated by the electronic 
structure. This hypothesis, as mentioned in the Intro- 
duction, is confirmed by the shell structure adopted by 
the atoms in thin nanowires-^ One of the main draw- 
backs of this assumption is that it neglects any possible 
back action of the atomic structure on the electronic mo- 
tion. However, these effects, though hard to evaluate, are 
expected to be minimal for the wires considered in this 
paper, which are thin enough that crystalline structure is 
not expected — and is indeed not observedi*^ — but thick 
enough for a continuum approximation to be reasonable. 

The simplicity of the model makes simulations of large 
systems over long times tractable. The initial wire for 
the simulation of Sec. |TT] corresponds to 4 • 10'' atoms, 
over an evolution time of seconds. This is far beyond 
what any simulation based on more "realistic" mod- 
els, such as ab initio method o^^'^" or classical molecu- 
lar dynamics, ^'^^ can hope to achieve in the foreseeable 
future, as they are currently limited to at most a few 
nanoseconds. Furthermore, classical molecular-dynamics 
simulations, which can treat relatively large systems, do 
not include electron-shell effects and are therefore unable 
to even stabilize a long nanocylinder. The NFEM — with 
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all its simplifying assumptions, and thanks to them — 
is to date the only model allowing realistic predictions 
of the long-time dynamics of metal nanowires. As the 
latter is now accessible experimentally with atomic reso- 
lution, these predictions are falsifiable and can be rigor- 
ously tested. 

Finally, although the NFEM is, in principle, limited to 
simple free-electron-like metals, such as the alkali met- 
als and to some extent noble metals, confinement of the 
electronic transverse motion is a very general feature of 
metals as soon as the transverse linear dimension of the 
system is of the order of the Fermi wavelength. As such. 



its consequences, from the magic radii for the stability to 
the discreteness of the ionic dynamics, are expected to 
be quite general. 
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